ACEEE Int. J. on Electrical and Power Engineering, Vol. 01, No. 03, Dec 2010 



Biogeography Based Optimization Approach for 
Optimal Power Flow Problem Considering Valve 

Loading Effects 

P. K.Roy \ S. P. Ghoshal 2 , S.S.Thakur 2 
1 Dr. B. C. Roy Engineering College, Department of Electrical Engineering ,Durgapur, West Bengal, India 

Email: roy_provas@yahoo.com 

2 National Institute of Technolog, Department of Electrical Engineering, Durgapur, West Bengal, India 

Email: {spghoshalnitdgp@gmail.com, sst_nit_ee@yahoo.co.in} 



Abstract — This paper presents a novel Biogeography Based 
Optimization (BBO) algorithm for solving multi- objective 
constrained optimal power flow problems in power system. In 
this paper, the feasibility of the proposed algorithm is 
demonstrated for IEEE 30-bus system with three different 
objective functions and it is compared to other well 
established population based optimization techniques. A 
comparison of simulation results reveals better solution 
quality and computation efficiency of the proposed algorithm 
over particle swarm optimization (PSO), Real Coded Genetic 
algorithm (RGA) for the global optimization of multi- 
objective constrained OPF problems. 

Index Terms — Biogeography, Optimal Power Flow, Particle 
Swarm Optimization, Genetic Algorithm, Evolutionary 
Programming, Mutation, Migration. 

I. Introduction 

Optimal Power Flow (OPF) problem may be defined as 
determining the optimal settings of a given power system 
network that optimizes a certain objective function while 
satisfying power flow equations and inequality constraints 
like system security and equipments operating limits. 
Several conventional optimization techniques such as linear 
programming, interior point method, reduced gradient 
method and Newton method [1, 2], have been applied to 
solve OPF problem assuming continuous, differentiable 
and monotonically increasing cost function. However, 
higher order nonlinearities and discontinuities are observed 
in practical input-output characteristics due to valve point 
loading. Conventional methods have failed in handling 
nonconvex and nondifferentiable optimization problems. In 
the recent years, many evolutionary algorithms (EA) like 
GA [3-5] simulated annealing (SA) [6], particle swarm 
optimization (PSO) [7-10], EP [11], and hybrid 
evolutionary programming (HEP) [12], ant colony 
optimization (ACO) [13] and Bacteria foraging 
optimization (BFO) [ 14] have been proposed to solve non- 
convex, non-smooth and non-differentiable optimization 
problems. 

Amongst the above population based algorithms, the 
annealing schedule of SA should be tuned carefully; 
otherwise it may produce suboptimal solutions. The GA 
method is usually faster than SA method because GA has 
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parallel search technique. Traditional GA also differs from 
EP in two aspects; EP primarily relies on mutation and 
selection, but no crossover like traditional GA and EP uses 
the real values of control parameters but not their coding as 
in traditional GA. Hence, considerable computation time 
may be saved in EP. Real coded GA (RGA) [5] has been 
introduced to solve the OPF problems more efficiently with 
significant reduction in the computation time. ACO is 
based on foraging behavior of ant species. Solution 
candidates, called ants in ACO, communicate with other 
members of the ant colony by depositing pheromone to 
mark a path. High concentrations of pheromones indicate 
more favorable paths that other members should follow in 
order to reach the optimal solution. BFO is a bio-inspired 
technique, applied to solve power system optimization 
problems by Ghoshal et al. [14], but its optimization time is 
very high. 

Biogeography based optimization (BBO) [15] is a new 
optimization algorithm and it has never been used in power 
system optimization so far. The application of 
biogeography to optimization has been first presented in 
[15] and it describes how a natural process can be modelled 
to solve general optimization problems. 

This paper presents BBO algorithm for solving multi- 
objective OPF problem of IEEE 30-bus system. The 
simulation results of BBO algorithm are compared to the 
results of PSO, and RGA and the computation efficiency is 
found to be more efficient in BBO than other conventional 
methods. 



II. Mathematical problem formulation 



A 



Objective Function 

In this paper, three different objective functions are 
considered to determine the effectiveness of the proposed 
algorithm. The objective functions are as follow: 

1) Minimization of Fuel Cost 

The aim of this objective function is to minimize the 
total fuel cost having valve point effect, while satisfying all 
equality and inequality constraints and it is formulated as: 
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(i) 

2) Minimization of Fuel Cost and voltage deviation 

The objective of this type of problem is to minimize the 
voltage deviations of buses and fuel cost of the generating 
units and may be mathematically modelled as follows [2]: 
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3) Minimization of Fuel cost and Transmission Loss 

To reduce the transmission losses in addition to fuel cost, 

mathematical formulation of this objective function is 

formulated by (3) [2]: 
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Subject to the following equality and inequality 
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where Cl i , t> ( - , C ( - , U i , ^ ( - : Cost coefficients of ith unit, 

■*pj » i^g/ » -*z,j 5 \LlI '■ Active & reactive power 
generation and load of i-th unit, 
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:Minimum and maximum generator voltages, load voltage, 
active and reactive power generations, tap setting and 
VAR injection of i-th unit respectively, 



line, 



P loss , P d , P di : Total loss, total load and load of i-th bus, 
NB,NG,NL,NC,NTL,NT : Number of buses, 
generator buses, load buses, shunt compensators, 
transmission lines and tap setting transformers. 

III. Algorithms 

A. Particle swarm optimization 

In a PSO system, multiple candidate solutions coexist 
and collaborate simultaneously. Here each candidate 
solution is associated with a velocity. Each candidate 
solution, called a particle, flies in the problem space 
(similar to the search process for food of a bird swarm) 
looking for the optimal position. A particle with time 
adjusts its position according to its own memory and 
cognitive experience, while adjusting to the social 
experience of neighbouring particles as well. If a particle 
discovers a promising new solution, all the other particles 
will move closer to it, exploring the region more 
thoroughly in the process requiring less computational 
book keeping and generally few lines of code. Based on 
PSO concept, mathematical equations for the searching 
process are: 

Velocity updating equation: 

v* +1 =wxvf +clxrlx(pBes?-Xi)+c2xr2x(gBest-Xi) (6) 



£+1 k k+\ 

Position updating equation: X t = X t + V t 



(V) 



where x, : current position of the i-th particle at k-th 
iteration, x t : modified position of the of i-th particle, V i 

k+\ 

current velocity of the i-th particle at k-th iteration, V i 
modified velocity of the of i-th particle, cl,c2 
weighting factors, r\, r2 : random numbers in the interval 
[0, 1], 

pBestj : position best of i-th particle, gBest : group 

best, W : inertia weight factor. 

Main steps of PSO techniques are given below: 

• Initialization of floating point particles of n„ 

population, each consisting of the variable 
parameters, 

• Evaluation of objective function value for each 
particle, 

• Search for individual minimum objective 
function value and corresponding individual 
best particle, 

• Search for global minimum objective function 
value and its corresponding global best particle, 

• Selection of elite particles like GA's selection 
or Taguchi selection, copying them over the 
unselected particles, 
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Velocity updating, 

Position updating, 

Individual best particle updating, 

Global best particle updating, 

Iteration updating and stopping criteria. 

B. Real Coded Genetic algorithm (RGA) 

Traditional binary coded GA suffers from few 
drawbacks when applied to multi-dimensional and high- 
precision numerical problems. The situation can be 
improved if GA is used with real number data. Different 
steps of RGA are 

Real coded initialization of each chromosome. 

Selection operation based on computation of 

fitnesses and merit ordering 

Crossover and Mutation operation 

Sorting of the fitness values in increasing order 

among parents and off-springs 

Selection of the better chromosomes as parents 

of the next generation 

Updating of genetic cycle and stopping 

criterion. 

C. Biogeography Based Optimization (BBO) 

1) Overview of BBO technique 

BBO [15] has been developed based on the theory of 
Biogeography [16]. BBO concept is mainly based on 
Migration and Mutation. The concept and mathematical 
formulation of Migration and Mutation steps are given 
below: 

Migration: BBO algorithm [15] is similar to other 
population based optimization techniques where population 
of candidate solutions can be represented as vectors of real 
numbers. Each real number in the vector is considered as 
one suitability index variable (SIV). Fitness (in BBO, a 
term called habitat suitability index (HSI)) of each 
candidate solution is evaluated using its SIVs. HSI 
represents the quality of each candidate solution. High HSI 
solution represents better quality solution and low HSI 
solution represents inferior solution in the optimization 
problem. The emigration and immigration rates of each 
solution are used to probabilistically share information 

between habitats. Immigration rate (A s ) and emigration 
rate ( jil s ) can be evaluated by (8) and (9) [15]: 
f \ 



A 



1 
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(9) 



where / : maximum immigration rate, 
£ : maximum emigration rate, 
S : nunber of species, 

S max : maximum number of species, 
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Using Habitat Modification Probability each solution is 

modified based on other solutions. Immigration rate, X s of 

each solution is used to probabilistically decide whether or 
not to modify each SIV in that solution. If a SIV in a given 

solution is selected for modification, emigration rates, jU^ 

of other solutions are used to probabilistically select which 
of the solutions should migrate a randomly selected SIV to 
that solution. The main difference between recombination 
approach of evolutionary strategies (ES) and migration 
process of BBO is that in ES, global recombination process 
is used to create a completely new solution, while in BBO, 
migration is used to bring changes within the existing 
solutions. In order to prevent the best solutions from being 
corrupted by the immigration process, few elite solutions 
are kept in BBO algorithm. 

Mutation: Due to some natural calamities or other events 
HSI of a natural habitat can change suddenly and it may 
deviate from its equilibrium value. In BBO, this event is 
represented by the mutation of SIV. Species count 
probabilities are used to determine mutation rates. The 
species count probability can be calculated using (10) [15]. 



-i\+Ms)Ps+Ms + l's +1 



- iK + ft )Ps + Vl*J-l + ft + A + 1 l ^ S ^ S,™ " 1 

-a+ft)^+Vi^-i s = s ma 

(10) 
Each population member has an associated probability, 
which indicates the likelihood that it exists as a solution for 
a given problem. If the probability of a given solution is 
very low then that solution is likely to mutate to some other 
solution. Similarly if the probability of some other solution 
is high then that solution has very little chance to mutate. 
Therefore, very high HSI solutions and very low HSI 
solutions are equally improbable for mutation i.e. they have 
less chances to produce more improved SIVs in the later 
stage. But medium HSI solutions have better chances to 
create much better solutions after mutation operation. 
Mutation rate of each set of solution can be calculated in 
terms of species count probability using the equation (11) 
[15]: 



(')■ 



r l~P s " 
P 

\ max / 



(ii) 



where m(s) :the mutation rate for habitat having S 

species, 
(8) 

m max • Maximum mutation rate, 

P max : Maximum probability. 

This mutation scheme tends to increase diversity among 
the populations. Without this modification, the highly 
probable solutions will tend to be more dominant in the 
population. This mutation approach makes both low and 
high HSI solutions likely to mutate, which gives a chance 
of improving both types of solutions in comparison to their 
earlier values. Few elite solutions are kept in mutation 
process to save the features of a solution, so if a solution 
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becomes inferior after mutation process then previous 
solution (solution of that set before mutation) can revert 
back to that place again if needed. So, mutation operation is 
a high-risk process. It is normally applied to both poor and 
good solutions. Since medium quality solutions are in 
improving stage so it is better not to apply mutation on 
medium quality solutions. 

Here, mutation of a selected solution is performed 
simply by replacing it with randomly generated new 
solution set. Other than this any other mutation scheme that 
has been implemented for GA can also be implemented for 
BBO. 

2) BBO algorithm applied to OPF 

In BBO algorithm, applied to the OPF problem, each 
candidate solution (i.e. each habitat) is defined by a vector 

with m SIVs where each SIV represents the value of 
independent variables such as active powers of all 
generators except slack bus, generators' voltages, tap 
settings of regulating transformers and reactive power 
injections. If there are m independent variables, then the i- 
th individual habitat H' can be defined as follows: 

H' = SIV '" = [SIV " , SIV n ,...., SIV "" ] (12) 

where i — \,2,...,n;q — \,2,...m 

The algorithm of this method is as enumerated below. 

Step 1 : Initialization of the BBO parameters. 

Step 2; The initial positions of SIVs (i.e. independent 
variables such as active powers of all generators except 
slack bus, generators' voltages, tap settings of regulating 
transformers and reactive power injections) of each habitat 
should be randomly selected while satisfying different 
equality and inequality constraints of OPF problems. 
Several numbers of habitats depending upon the population 
size are being generated. Each habitat represents a potential 
solution. It is similar to a chromosome string in GA. 

Step 3: Perform load flow using any classical technique 
and determine all dependent variables such as load 
voltages, active power of slack bus, generators' reactive 
powers etc. In this paper load flow is performed by 
Newton-Raphson method. 

Step 4: Calculate each HSf i.e. value of objective 
function for each i-th habitat of the population set n for 

given emigration rate jU s , immigration rate X s and 

species S . 

Step 5: Based on the HSI values elite habitats are 
identified. 

Step 6: Each non-elite habitat is modified by 
performing probabilistically immigration and emigration 
operation as described in section III. C.l. 

Step 7: Species count probability of each habitat is 
updated using (10). Mutation operation is performed on the 
non-elite habitat and HSI value of each new habitat is 
computed. 

Step 8: Feasibility of a problem solution is verified i.e. 
each SIV should satisfy equality and inequality constraints. 

Step 9: Go to step 3 for the next iteration. 
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Step 10: Stop iteration. 

IV. Input Parameters 

After several runs, the following input control 
parameters are found to be best for optimal performance of 
the proposed algorithm. Habitat Modification Probability, 
p mod =1; Mutation Probability = 0.005, maximum 
immigration rate,/ =1, maximum emigration rate, £ =1, 
step size for numerical integration, dt =1, maximum 
iteration cycles = 100, elitism parameter = 4, number of 
SIVs = number of generator units and number of Habitats 
= population size =50. 

V. Simulation results and discussions 

The applicability and validity of all the algorithms have 
been tested on IEEE 30-bus system for multi-objective 
OPF problems. This system contains 6 generating units 
connected in bus no. 1, 2, 5, 8, 11 & 13; 30 buses, 4 
regulating transformers connected between the line 
numbers 6-9, 6-10, 4-12 & 27-28; 2 shunt compensators 
connected in bus no. 10 & 14; 41 transmission lines. The 
bus data, transmission line data and load data of this system 
are taken from [17]. Generating unit capacity and cost 
coefficients data gmmd$]^$d&kfy-Wlmmpmtmtnfi&iable of the 
of all buses are considered within the range of 
[0.95/?.m,1.05/?.m.]- Tap settings of regulating 
transformers are within the range of [0.9 p.u., 1 .1 p.u.]- 
The VAR injection of the shunt capacitors are taken within 
the interval of [0 p.u., 0.3 p.u.]. The programming has been 
written in MATLAB language, using MATLAB 7.1 on 
core 2 duo processor, 2.0 GHz with 1GB RAM. 

The simulation results and convergence characteristics 
with three different objective functions namely, 
minimization of fuel cost, minimization of voltage 
deviation together with fuel cost and minimization of cost 
as well as transmission loss are given in Table II, Table III, 
Table IV, Fig. 1, Fig. 2 and Fig. 3 respectively. The active 
power, reactive power, bus voltage, VAR injection, tap 
settings, transmission loss (TL), summation of voltage 
deviation ( S vD = V |v.-l|) are on p.u. basis, fuel costs 

(FT) are in $/hr and computational times per iteration (CT) 
are in seconds . The simulation results show that BBO 
outperforms PSO and RGA in terms of solution quality 
(objective function, OF) and computation efficiency 
without violating any operating constraints. 

Table I 
genarating unit capacity and cost coefficients data 



Bus 
No. 


P . 

x mm 

p.u. 


p 

* max 

p.u. 


vtmin 

p.u. 


L/max a 

p.u. $/h 


b c d 

S/MW $/MW 2 h$/h 
h 


e 
rad/MW 


1 

2 







0.8 

0.8 




-0.4 


0.1 0.02 
0.5 0.0175 


2 
1.75 







300 

200 


0.2 
0.22 


5 





0.4 


-0.4 


0.4 0.0250 


3 





200 


0.35 


8 





0.5 


-0.1 


0.4 0.0625 


1 





150 


0.42 


11 





0.3 


-0.06 


0.24 0.025 


3 





200 


0.35 


13 





0.55 


-0.06 


0.24 0.00834 


3.25 





100 


0.3 
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Table II 
Simulation results of case 1 with different algorithms 



Alg. 



PSO 



RGA 



BBO Alg. PSO RGA 



BBO 



PG, 
PG: 
PG 5 
PG 8 
PG„ 
PG,: 
Vi 

v 2 
v 5 
v 8 

Vn 

V13 



0.7819 
0.5707 
0.3569 
0.3782 
0.2689 
0.5277 
1.0019 
0.9849 
0.9678 
0.9932 
1.0414 
1.0789 



0.7845 
0.5691 
0.3601 
0.3811 
0.2696 
0.5216 
0.9923 
0.9761 
0.9671 
0.9912 
1.0328 
1.0381 



0.6283 
0.7153 
0.3592 
0.4500 
0.1795 
0.5408 
1.0471 
1.0365 
1.0165 
1.0285 
1.0931 
1.0974 



QCio 
QC 24 

TC M 
TC(,_io 
TC 4 -, 2 
TC27-28 

OF 

FC 

TL 

SVD 

CT 



0.2964 
0.0839 
1.0316 
0.9329 
0.9275 
1.0373 
1074.5 
1074.5 
0.0503 
0.8073 
14.201 



0.2764 
0.1038 
1.0661 
0.9387 
0.9851 
1.0683 
1066.1 
1066.1 
0.0520 
0.9010 
12.653 



0.2894 
0.1012 
0.9600 
0.9652 
0.9235 
0.9501 
1058.7 
1058.7 
0.0394 
1.7648 
13.391 




RGA 
-PSO 

- mo 



10 20 30 40 50 60 70 80 90 100 

No of Iterations 



Fig. 1 Convergence Characteristics with different Algorithms (Case 1) 

Table III 
Simulation results of case 1 with different algorithms 



Alg. 

PG, 
PG : 
PG, 
PG 8 

PG,, 



PSO 

0.6247 
0.7111 
0.2741 
0.4619 



RGA 

0.6213 
0.7186 
0.2710 
0.4597 



BBO 

0.6253 
0.7211 
0.2695 
0.4687 



Alg. 

QC„, 
QC 24 
TC M 
TC 6 - 



PSO 

0.0911 
0.1915 
0.9588 
1.0329 



RGA 

0.0010 
0.1729 
1.0153 
0.9000 



BBO 

0.2921 
0.2863 
0.9974 

1.0274 



0.2690 0.2709 0.2683 TC 4 . 1.0065 1.0055 0.9873 



PG, 3 


0.5402 
1.0110 
0.9989 


0.5389 
1.0150 
1.0011 


0.5295 
1.0119 
1.0013 


TC 27 - 21 . 


0.9473 


0.9468 


0.9668 


V, 


OF 


1200.9 


1197.4 


1192.6 


V 2 


FC 


1191.3 


1190.4 


1184.1 


v 5 


0.9918 


0.9941 


0.9918 


TL 


0.0470 


0.0464 


0.0484 


V 8 


0.9999 


1.0010 


1.0013 


SVD 


0.3273 


0.2582 


0.2333 


V„ 


1.0401 


1.0334 


1.0123 


CT 


14.237 


12.785 


13.573 


V, 3 


1.0155 


1.0391 


1.0142 












No of Iterations 

2 Convergence Characteristics with different Algorithms (Case 2) 

Table IV 
Simulation results of case 3 with different algorithms 



Alg. 



PSO 



BBO 



Alg. 



PSO 



RGA 



BBO 



RGA 



PG, 
PG: 

PG 5 

PG 8 

PG,, 
PG, 3 

V, 

v 2 

Vs 
V 8 
Vn 

V,3 



0.6244 
0.7221 

0.3614 

0.3668 

0.2661 

0.5313 

1.0559 

1.0480 

1.0266 

1.0409 

1.0986 

1.0911 



0.7703 
0.5779 
0.3604 
0.3711 
0.2695 
0.5231 
.1000 
1.0883 
1.0836 
1.0901 
1.1000 
1.1000 



0.6341 
0.7005 

0.3599 

0.3707 

0.2698 

0.5347 

1.1000 

1.0905 

1.0705 

1.0851 

1.1000 

1.1000 



QCio 

QC 24 

TC M 

TCe-io 

TC 4 _, 2 

TC27-28 

OF 

FC 
TL 
SVD 
CT 



0.2938 
0.1096 

0.9130 

1.0375 

0.9204 

0.9237 

1156.7 

1141.3 

0.0381 

2.1289 

14.318 



0.3000 
0.1045 

0.9782 

1.1000 

0.9681 

0.9749 

1152.5 

1137.1 

0.0383 

2.4912 

12.814 



0.3000 
0.0995 

1.1000 

0.9202 

0.9541 

0.9692 

1146.6 

1132.3 

0.0357 

2.5316 

13.579 



^ 




RGA 

PSO 

BBO - 







No ol Iterations 

Fig. 3 Convergence Characteristics with different Algorithms (Case 3) 

Conclusions 

In this paper a novel Biogeography Based Optimization 
algorithm is proposed to solve multi-objective OPF 
problems. The test results clearly demonstrate that BBO 
outperforms other reported methods in terms of solution 
quality, computational efficiency, convergence 
characteristics for both IEEE 30-bus system. It is revealed 
that among all the algorithms, Biogeography Based 
Optimization is capable of achieving global optimal 
solution. 
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